Uniqueness of the equation for state-vector collapse 
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The linearity of quantum mechanics leads, under the assumption that the wave function of- 
fers a complete description of reality, to grotesque situations famously known as Schrodinger's 
cat. Ways out are either adding elements of reality or replacing the linear evolution by a 
nonlinear one. Models of spontaneous wave function collapses took the latter path. The 
way such models are constructed leaves the question, whether such models are in some sense 
unique, i.e. whether the nonlinear equations replacing Schrodinger's equation, are uniquely 
determined as collapse equations. First Adler, and more recently Weinberg, worked on iden- 
tifying the class of nonlinear modifications of the Schrodinger equation, compatible with 
general physical requirements. Here we argue that such considerations lead in fact to the 
usual collapse models. 
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INTRODUCTION 



The superposition principle of quantum mechanics has always been an open issue. The strongest 
reason is that the collapse of the wave function, which has been formulated to avoid macroscopic 
superpositions, is clearly an artifact |li]. To overcome this situation, collapse models have been 
devised, which incorporate wave-function's collapse, together with the Schrodinger evolution, into 
one single (nonlinear and stochastic) dynamical law {2I-8]. Independently from this, Weinberg ques- 
tioned whether quantum mechanics, being a linear theory, should be regarded as an approximation 
of an underlying non-linear theory Q], pretty much in the same way in which Newtonian gravity 
is a weak-limit limit of general relativity. Moreover when one combines quantum mechanics and 
gravity even in the most naive way, by taking into account the gravitational interaction of the wave 
function with itself, one obtains nonlinear evolutions 1CH121] . 



Quantum linearity has been challenged also experimentally, and the interest and efforts are 
constantly increasing. Perhaps the most famous tests are diffraction experiments with macro- 
molecules 



13l4l5||. In the most recent experiment of this kind, quantum coherence was proven to 
hold for molecules with up to about 7000 amu [15]. The question remains open whether larger 
systems still enjoy this property. Other proposals, representing a direct test of quantum linearity, 
include satellite atom (Cs) [lfj], micro-mirrors nano-sphere interferometric experiments [l8|. 
To these, one should add the vast zoo of indirect tests of the superposition principle {(], [ll]]. 

Among all proposals for nonlinear modifications of quantum mechanics, only collapse models 
survived. The most important reason is that, in response to Weinberg's attempt to modifying 



quantum theory in a nonlinear sense, N. Gisin 



20 



211 ] . followed by J. Polcinski 22], proved that 



any nonlinear and deterministic modification of the theory yields superluminal signalling. Therefore 
only stochastic nonlinear modifications are allowed. Collapse models fall in this class. However, it 
is not clear at all whether they should be the only possible models, since their nonlinear structure 
is very specific and given by the following stochastic differential equation ij: 



iHdt + Y,(L k - £ ktt ) dW ktt ( L k L k ~ 2 ^,tL k + \l k , t \ 2 ) dt 



fc=l 



k=l 



4>t, 



4 jt = I(^,(4 + J L fe )^} 



(i) 



(2) 



where H is the standard quantum Hamiltonian of the system, L k are linear operators defining the 
preferred basis for the collapse, and W k % standard Wiener processes. As one can see, nonlinearity 
enters in a very specific way. Of course the above equation can be generalized, e.g. by by considering 



other kinds of noise. One can also consider jump-type models, like the GRW model Q], instead of 
diffusion processes. Nevertheless, the basic structure does not change. 

S.L. Adler and T.A. Brun were the first to raise the question whether indeed collapse models are 
the only models compatible with the no-faster-than-light requirement 
used in the context of Gisin's theorem 
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241 ] . This assumption is 



21] and implies that, whatever the evolution for the state- 



251 ] . Such considerations 



261 ] and of 



vector, the equivalence among statistical ensembles must be preserved 
lead to linear evolutions of statistical operators. According to the work of G. Lindblad 
V. Gorini, A. Kossakowsky and E.C.G. Sudarshan [^tJ, they have to be of the Lindblad type The 
question can be thus phrased as follows: 
Question: Given a diffusion process for the wave-function: 



#t = A{^ t )dt + J2 B k (rp t )dW Kl 



(3) 



k=l 



where A(ip) and B k {ip) are unspecified nonlinear operators, find the conditions under which it 
generates a Lindblad type of equation: 



dpt sr*- ( t 1 t 1 + \ 

— = -i[H,pt\ + 2^ \ L kPtL k - -L^Lkpt - -p t L' k L k I 

k=l ^ " - / 



(4) 



for the statistical operator pt = M[\i()t){ipt\]. 

One should specify that the above quoted theorems of Lindblad and of Gorini, Kossakowski and 
Sudarshan require other assumptions besides linearity of the dynamics; in particular, they require 
a quantum-dynamical semi-group type of evolution, and complete positivity. The first assumption 
basically amounts to Markovianity of the dynamics which usually is accepted for fundamental 
evolutions as we consider here. The second assumption, can be motivated by the requirement that 
reduction to subsystem descriptions, i.e. tracing out of "bath" variables should produce equations 
of the same type as the initial one. With the above specifications, one can say that the no-faster- 
than-light condition implies a Lindblad type of evolution for the statistical operator. 

Adler and Brun provided a partial answer to the question, their proof being restricted to the 
very special case of linear operators 28j]. Recently, also S. Weinberg [8| raised the same question. 



Here we answer the question in full generality, by proving the following statement. 

Main Result. Eq. ([3]) leads to Eq. ([4]) if and only if the operators A(tjj) and B k {^ijj) take the 

following form: 

1 n 

AW = ~iH-W (L{L k - 2$4*> + \tf} | 2 ) (5) 



B k ($) = L 



k=l 



l k,t » 



l k,t 



1 



«0> 



(6) 
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modulo unimportant global phase factors. The operators L, and their "adjoints" are defined as 
follows, 

n n 

L? ] ■= E «*(V0 Lj, L™ : = £ 4" ( 7 ) 

where Ukj{ip) are the unknown coefficients of a n x n unitary matrix. This unitary freedom, which 
is a symmetry of the Lindblad equation © but not of the diffusion equation ([3]), implies that the 
allowed diffusions range from linear diffusion (no collapse) to the usual collapse equation. The 
collapse rate for the diffusion is maximal for the usual collapse equation ([I]) . 

This is a remarkable result: it basically says that collapse models are the only possible nonlinear 
extensions of the Schrodinger equation, compatible with the no-faster-than-light assumption. 



II. PROOF OF THE STATEMENT 



First of all let us note that, without loss of generality, in Eqs. (J3j> and (j3J) we can take m = n. 
In fact if, for example, m > n, then we can add m — n null operators to Eq. @ to increase the 
sum up to m. 

Step 1. By using Ito calculus [29], one can show that Eq. (|3|) leads to the following master equation 
for p t :=E[|Vt><^|]: 

^ = E (\A(^t))(A\ + \^t)(A^ t )\ + p \B k (i> t ))(B k & t )\\ . (8) 

Let us choose a deterministic initial value ip in Eq. §3$). Correspondingly, we set po := |V ; )(V ; I as 
the initial value in Eq. (J3J). Let us set both Eq. (JH]) and (jl]) equal at time t = 0. This will be 
sufficient, because the equations are Markovian and time-homogenous. We have: 



\AM)ei>\ + \i>)(Aw\ + E \B k (4>))(B k (n = 

k=l 

n / i i \ 

= | - iH^M + m-%H^\ + E - -\L\L^)^\ - -1^(4^1 • (9) 



k=l v y 

Step 2. Since B k and L k are the only operators appearing in the "bra" and "ket" part simulta- 
neously, then the B k are essentially fixed by the L k . To find out how, we restrict Eq. to the 
orthogonal complement of tp. Writing the corresponding orthogonal decomposition of vectors as 
v = v± + U|| , we get: 

n n 

E \B k W±)(B k &)±\ = E . (10) 

k=l k=l 
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This implies that: 

n 

B k (i/>)± = Y, u kjW(Lji>)±, (11) 

J'=l 

where the functionals u k j(ip) are the coefficients of a '(/'-dependent unitary matrix [3^]. Therefore 
we can write: 

n 

B k w = y, u kM ( L ^)± - (12) 

for unspecified functionals 1^, of V>- Equivalently, we can write: 

n 

B k {^) = Y ukM)h1> - 4^ = - efi> (13) 



where LjJ is defined as ([7|), and we have defined: 



S^ep 5. We now determine the operator A. For this, we substitute the previous result into Eq. Q 
and we multiply with ip from the left and v -Lip from the right. We get: 

n 1 n 

(Affl,v) - Y^i^M = i(H^v) --YiLlL^v). (15) 
k=l fc=i 

Since u is a generic vector orthogonal to ip, the above equation implies that the parts of A(ip) — 

Ylk=i ^*L^rp and of —iH — \ Y^k=i ^ k ^ J k' l P orthogonal to ip coincide, while the parallel part of 

A(ip) is arbitrary. Therefore, similarly to the previous step, one can write: 

1 - 

A{iP) = -zHiP - i £ (l) k L^ - 2£f*L^V) + Gty)* (16) 

k=l 

with arbitrary complex functional G. 

Step 4- Substituting Eqs. (fT3j) and ([TBI) into ([9]) and multiplying from both sides with ip gives: 



G^) = -\Y.\^ ] \ 2 + ^) (17) 

k 

with real functional g. 

Step 5. Now we apply the condition of norm preservation for the diffusion process given by Eq. ([3]) . 
The stochastic differential equation for the square norm of ipt is: 

(n \ n 

(iMM) + (A(V0,V) +^(S fc (^), J B fc (V)) dt + J3«^B fc W) + (B k ty),$))W k>t 
k=i J k=i 

(18) 
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which should be zero, for the norm to be conserved. Because of the previous relations, one can 
show that drift term is automatically zero, while requiring that also the diffusion terms are zero, 
give: 

(^4^) + <4^V,^}-2Re[4^] =0, k = l,...n (19) 

or: 

^) = I^ 5 (4^ +L w w _-^ W) (20) 

for arbitrary real- valued functionals hf~(ip). Eq. (fT3j) , with £^ given by Eq. ([20]) . identifies the 
operators B^^i/j). Similarly, Eq. (fT6|) . with G(tp) given by Eq. (fTT|) . identifies the operator A(ip). 
So far, the unspecified quantities are the functionals h^ip) and g(ip), as well as the coefficients 
u k j(ip). 

It is easy to see that he functionals g and hk are unimportant global phase factors. In fact, 
tpt := e^^o9( s ) ds ^ t satisfies Eq. (|8|) with g = 0. In a similar way, <j> t := e~*^=i fo h k(s)dW k , a ^ 
satisfies Eq. ([5]) with g = = 0, as a simple application of Ito's formula shows. Therefore the 
only degree of freedom is encoded in the the unitary matrix with coefficients Ukj(ijj)- 

By straightforward use of Ito calculus, the operators A and B^ we have found, indeed reproduce 
the Lindblad equation (j4j). Therefore also the 'if part is proven. This concludes the proof of the 
statement. 

To understand the meaning of the unitary freedom, let us consider the simpler case of only one 
self-adjoint Lindblad operator L. In such a case, the unitary matrix reduces to a complex phase 
factor e^M. Eq. ©, with the previous specifications for A and B^, becomes: 

dip t = -iH^ t dt+{e lf ^L^ t - (^,L^)^cos(/(Vt))) dW t 

~ (L 2 lH - 2e^^)(^,L^)L^cos(/(Vt)) + \(ipt, i^)| 2 ^cos 2 (/(^))) dt, (21) 

When / = 0, one recovers the standard collapse equation [5j. For / = ir/2 instead, all non- 
linear terms vanish: one obtains a linear and random evolution, which, written in the Stratonovich 
formalism — which is closer to the physicists' usual formalism — becomes: 

= -i[H - Lw t ]iP t , (22) 

with Wt = dWj/dt. This is a standard Schrodinger with a random potential, which does not give 
rise to the collapse of the wave function (since it is linear), but nevertheless originates the same 
Lindblad equation for the density matrix. This is an entirely expected result, and the previous 
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literature is aware of such a fact, i.e. that there are infinite stochastic unravellings of the same 
master equation 



31 



32] 



One way to understand the meaning of / is to say that it measures the "anti-hermitian" part of 
the coupling of the noise with the wave function, which is the one responsible for the collapse. The 
maximum such coupling is obtained precisely for / = 0. It is in fact, let us neglect the Hamiltonian 
H and let us compute the equation for the variance Vt '■= (ipt, L 2 ipt) — (ipt, Lipt) 2 , which is a measure 
of how the stavector is delocalized over the eigenstates of L. One obtains [33]: 

dV t = -4cos 2 (/(^M 2 dt + 2cos(/(^))(Vt, (L - fa, L^ t )fi> t )&W t . (23) 

The drift term (which proportional to — V 2 ) drives Vt to zero for large times and for almost 
any realization of the noise, as expected from a collapse equation. The rate is proportional to 
cos 2 (/(V't)) and of course is maximal for / = 0. This is the standard choice for collapse models. 



III. FOKKER-PLANCK EQUATION 



It is interesting to understand in which way the probabilities for the wave functions develop on 
Hilbert space which then allow to classify the possible diffusions in more detail. For example we 
note that some diffusions have the same diffusive term in the Fokker Planck equation governing the 
probability distribution on Hilbert space: The Fokker-Planck equation for the probability density p 
with respect to the "Hilbert space volume element" VvpnDip] (R and I denoting real and imaginary 
part) reads 



dt 



- E 

rne{R,I} 

+ 1 E 



a 



d 2 



dip(x)idip(y) 
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(x)dx 

(£>(V)pW0) 



l,m&{R,I} ' 

if the diffusion matrix D is defined via its real entries 



{x, y\ m dxdy , 



[D^)](x,y\ m = Y,lBk(Mx)i[B k m(y)m ■ 



fc=l 



More precisely, choosing an orthonormal basis of the Hilbert space and introducing the correspond- 
ing complex coordinates x% = (xi t R, x\j),X2, ■ ■ ■ , one can write 

dpt(x) sr^ , r . , , NS 1 s d 2 



Of 



E 

k=l,2,... 
m£{R,I} 



dx 



k,m 



{[A(x)] Km p{x)) + - ]T 



2 , . dxi m dx j n 

m,n£{R,I} 



[ D (x)])™p{x] 
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where the diffusion matrix D has entries 

Wljj = Y,l B k( X )krn[Bk(x)) jt n. 
fe=l 

Thus, within this calculus, one has to consider the Hilbert space as a real one with double dimension. 
In particular, the complex matrix \Bk){B] c \ can be computed from D, but not vice versa - the 
distribution of ipt is not fixed by Ylk \Bk){Bk\- 

Replacing the matrix u by ou, where o is real orthogonal, leaves not only \Bk)(Bk\, but also 
D invariant; in particular, for any real orthogonal u, the diffusion matrix is the same as for m j = <5y 
(the standard collapse equation). With complex entries, this need not be the case. The simplest 
example is the choice Uy = iSij which, assuming that the are self-adjoint, gives B^ip) = iL^ip 
and A(ip) = —iHip^ J2 k Lfjp, i.e. a linear equation reproducing Eq. fll]) and generalizing (|22|) . 

We remark that having introduced the probability distribution p allows us to write the density 
matrix more explicitly as the second moment of the probability distribution on the Hilbert space, 
namely 

Pt{x,y) = j ip(x)ip(y)T>[ip] . 
IV. BORN'S STATISTICAL LAW 



In the previous discussion, we never required that the collapse of the wave function should 
occur according to Born's statistical law, which is the standard rule of quantum mechanics. But 
indeed one can prove that, in measurement situations, collapse model predict that outcomes are 
distributed according to the Born rule 



24, 



34( . Accordingly, as first noted in 3, Q] , this rule is the 



only probability rule which is compatible with the requirements of no-faster-than-light signalling 
and norm preservation. 

Here we wish to comment on the reason why this is indeed the case. This can be understood by 
considering the following situation. Assume that a state vector tp suddenly collapses — because of a 
not otherwise specified collapse dynamics — to one of the eigenstates |n) of a (self-adjoint) operator 
L, with a probability distribution p n . The operator describing the collapse effect is: 



2^ Pn 1 1 PL/A 1 12 ' ^> 



with P n = \o n )(o n \. Such an operator is trace preserving (the analog of norm conservation). If we 
also require T to be a linear operator (as required by the no faster-than-light signalling condition), 
than the only possible choice is to take p n = | |-fn l^) 1 1 2 , which is the Born's statistical law. 
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